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Abstract 

Considering the magnetic moment and polarizabilitics of the nucleon we emphasize the need for 
relativistic chiral EFT calculations. Our relativistic calculations arc done via the forward-Compton- 
scattcring sum rules, thus ensuring the correct analytic properties. The results obtained in this 
way are equivalent to the usual loop calculations, provided no heavy-baryon expansion or any other 
manipulations which lead to a different analytic structure (e.g., infrared regularization) are made. 
The Baldin sum rule can directly be applied to calculate the sum of nucleon polarizabilitics. In 
contrast, the GDH sum rule is practically unsuitable for calculating the magnetic moments. The 
breakthrough is achieved by taking the derivatives of the sum rule with respect to the anomalous 
magnetic moment. As an example, we apply the derivative of the GDH sum rule to the calculation 
of the magnetic moment in QED and reproduce the famous Schwinger's correction from a tree-level 
cross-section calculation. As far as the nucleon properties are concerned, we focus on two issues: 
1) chiral behavior of the nucleon magnetic moment and 2) reconciliation of the chiral loop and 
A-resonance contributions to the nucleon magnetic polarizability. 
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1 GDH sum rule and its derivatives 



Consider the elastic scattering of a photon on a target with spin s (real Compton scattering). The 
forward-scattering amplitude of this process is characterized by 2s +1 scalar functions which depend on 
a single kinematic variable, e.g., the photon energy oj. In the low-energy limit each of these functions 
corresponds to an electromagnetic moment — charge, magnetic dipole, electric quadrupole, etc. - 
of the target. For example in the case of the nucleon, the forward Compton amplitude is generally 
written as 

T(u) = t ■ sf(u) + iS ■ (e* x g) g(u) , (1) 

where e, is the polarization vector of the incident and scattered photon, respectively; a are the 
Pauli matrices representing the dependence on the nucleon spin. The two scalar functions have the 
following low-energy expansion, 

/M = -^ M + (a E + p M )u J 2 + 0(u; i ) , (2a) 

<?M = -^^ + 7ou; 3 + 0(u, 5 ) , (2b) 

hence in the low-energy limit they are given in terms of the nucleon charge e and the anomalous 
magnetic moment (a.m.m.) k. The next-to- leading order terms are specified by the nucleon electric 
(ue), magnetic (Pm), and forward spin (70) polarizabilities. 

To derive the sum rules (SRs) for these quantities one assumes the scattering amplitude is an 
analytic function of u everywhere but the real axis 1 . This allows us to write the real parts of functions 
f{oj) and g{u)) as a dispersion integral of their imaginary parts. The latter, on the other hand, can be 
related to the total photoabsorption cross-sections by using the optical theorem, 

hn f(u) = — [<7i/ 2 (w) + ff3/a(w)] , (3a) 

Im g(u) = ^ [ay 2 {^) - 03/2H] , (3b) 

where o\ is the double-polarized total cross-section of the photoabsorption processes. Averaging over 
the polarization of initial particles gives the total unpolarized cross-section, o~t = \{o~\/2 + °3/2)- 

Finally one uses the crossing symmetry, meaning that the Compton amplitude of Eq. must be 
invariant under e' <-> e, 00 — > —u, and hence / is an even and g an odd function of energy:/(u;) = /(— oS), 
SM = -9{~u). 

Going through these steps one arrives at the following result (see, e.g., [1] for more details): 

1 f 00 °tW) . n 



/M = 5-5 / 72 2 ■ " ^' ■ ( 4 ) 

9^) = -in / -12 — 2 — ~ u duJ ' ( 5 ) 

with Act = (T 3 / 2 — o~i/2- These relations can be expanded in energy to obtain the SRs for the different 
static properties introduced in Eq. ([%]). In this way we can obtain the Baldin SR: 

ve + Pm = ^ f°^W, (6) 







1 Resonance poles may occur but lie on the second Riemann sheet. 
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the Gerasimov-Drell-Hearn (GDH) SR: 

e 2 K 2 ^ 1 roc A(J ^ 

2M 2 ~ n J ~ 
and a SR for the forward spin polarizability: 



duo , (7) 



1 f°° Aa(u) , 

As we all now know, impressive experimental programs to measure the total photoabsorption cross- 
sections of the nucleon have recently been carried out at ELS A and MAMI (for a review see Ref. [2]). 
These measurements are needed for an empirical test of the GDH SR, as well as for phenomenological 
estimates of og + @m and 70 via the other two SRs. The GDH SR is particularly interesting because 
both the left- and right-hand-side of this SR can reliably be measured, thus providing a test of the 
fundamental principles (such as unitarity and analyticity) which go into its derivation. 

Testing theories by using these SRs could also be fun and even useful in instances when the 
consistency of the theory is not transparent. Let us, for example, have a look at the left- and right- 
hand-sides of the GDH SR for the electron in QED. To lowest order in the fine-structure constant, 
a = e 2 /4-7r, the photoabsorption cross-section is given by the tree-level Compton scattering cross- 
section [3]: 



Acr(w) = — — 



2co 2 ( M\ , / 2uA 



+ 0(a 3 ). (9) 



On the other hand, the one-loop contribution to the electron a.m.m. is of order a and therefore the 
Ihs has no contribution of a 2 . Fortunately, the GDH integral over the tree- level cross-section Eq. Q 
vanishes, and thus, at this order, everything works out: 

= (10) 

as one could expect for such a fortunate theory as QED. At order a 3 the Ihs receives the contribution 
in the form of Schwinger's correction: n = a/2n. The calculation of cross sections at this order is 
quite a formidable task as it requires the knowledge of Compton scattering amplitude to one- loop, 
inclusion of the pair-production channel, and so on, cf. [4]. Instead, I want to follow a much simpler 
way to do calculations at this level [5]. 

Let us introduce a 'classical' (or 'trial') value of the a.m.m., kq. At the level of the Lagrangian 
this amounts to introducing a Pauli term for our spin-1/2 field: 

^Pauli = (*'Ko/4M) $ <V V , (11) 

where F is the electromagnetic field tensor and 07^ = (i/2)[y^,^ u ]. In the end we can put ko equal 
to zero, but for now the total value of the a.m.m. is k = kq + 5k, with 5k being the quantum effects. 
Note that 5k and the total cross-sections become explicitly dependent on kq. To get something new 
out of this we need to start taking derivatives of the GDH SR with respect to kq: 

roc j 

{An 2 a/M 2 )KK' = / A<t'(cj) — > (12) 
Jo w 

(4ir 2 a/M 2 )(K ,2 + KK") = / Aa"(uj) — , (13) 

Jo w 
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and so on. Now observe that to lowest order in a these relations simply read as 

(47r 2 a/M 2 )n^ n -^ = Aff W (w) — (14) 
Jo v 

where n is the order of the derivative with respect to Kq. This allows us in principle to compute k to 
order a k by using kth, (k — l)th, etc., derivatives of the cross-section computed, respectively, to order 
a k+1 , a k , etc. In this way, to lowest order we have the following sum rule: 

roc j 
(4vr 2 a/M 2 )^= / !_ — . (15) 



The striking feature of this sum rule is the linear relation between the a.m.m. and the photoabsorption 
cross section, in contrast to the GDH SR where the relation is quadratic. This restores the "balance of 
difficulty" in the two methods of calculating this quantity: the sum rule or the usual loop technique. 

Although the cross-section quantity Act' is not an observable, it is very clear how to determine it 
within a given theory. The first derivative of the tree-level cross-section with respect to Krj, at kq = 0, 
in QED takes the form [5]: 



a u si 2ira 2 
Aa'(io)\ n = --r- 

It is then not difficult to find that 



2Mu (_ 3M\ / 2uA 



6 -(mtw"V 2+ -J ln { 1+ M) 



(16) 



1 f \ i / \ i duo 2a 2 

- / Ad' (w) . = T7o- (17) 



Substituting this result in the linearized SR, Eq. (|15|) . we obtain k = a/2ir. Thus, the Schwinger's one- 
loop result is reproduced here by computing only a (derivative of the) tree-level Compton scattering 
cross-section and then performing a GDH integral. 



2 Magnetic moments and their chiral extrapolation 

Consider now the theory of nucleons interacting with pions via pseudovector coupling: 

9 



a 

/Lt" ) 



(18) 



where g is the pion-nucleon coupling constant, r a are isospin Pauli matrices, ir a is the isovector pion 
field. For our purposes this Lagrangian is sufficient to obtain the leading order results of chiral 
perturbation theory. 

To lowest order in g the photoabsorption cross section in this theory is dominated by the single 
pion photoproduction graphs as displayed in Fig.^ We find for the corresponding GDH cross sections: 



Aa {n+n) 
A ff (* - p) 



_kC_ 

M 2 x 2 
2ttC 
M 2 x 2 

0, 

2vrC 
M 2 x 2 



(2as + 1 - x) In 



a + A 
a — A 



2\[x(a - 2) + s(a + 2)] 



-li 2 In + 2X (s(3 ~ xa) 



(3- A 



-fj, in 



p-x 



+ (2as - 1 - x) In ^4 - 2s\ 
a — A 



(19a) 

(19b) 
(19c) 
(19d) 
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(a; i ^ H^L ^ h 



/ 



(b) —2^ Z_ 




Figure 1: Tree- level pion photoproduction graphs. The circled vertex corresponds to the Pauli cou- 
pling. 

where C = {eg/Air) 2 , fi = m n /M, is the pion mass, and 

s = M 2 + 2Mu>, s = s/M 2 , (20a) 
a = {s + M 2 - ml) /2s, (20b) 
P = (s - M 2 + m 2 n ) /2s = 1- a, (20c) 



A = (l/2s) V* - (M + m n ) 2 y/s — (M - m n ) 2 . (20d) 

As in the case of QED, the anomalous magnetic moment corrections start at 0(g 2 ), implying that 
the Ihs of the GDH SR begins at 0(g 4 ). Since the tree-level cross sections are 0(g 2 ), we must require 
that 

oo 

f — Aa {1) (u) = 0, for / = ir°p, vr+n, vr n, vr^p, (21) 
J u 

where Uth = m 7r(l +?ri 7r /2M) is the threshold of the pion photoproduction reaction. This requirement 
is indeed verified for the expressions given in Eq. — the consistency of GDH SR is maintained in 
this theory for each of the pion production channels. 

We now turn our attention to the linearized GDH sum rule. In this case we first introduce Pauli 
moments kq p and KQ n for the proton and the neutron, respectively. The dependence of the cross- 
sections on these quantities can generally be presented as: 

Aa(u; k 0p , K Qn ) = Aa{u) + K 0p Ao lp (u) + k „ Aa ln (u) 

+ Kl p Aa 2p (uj) + KQn A<7 2 nM + K 0p K 0n A(T lpln (oS) + . . . . (22) 

Furthermore, we introduce proton and neutron photoproduction cross sections A<jW and AcrW and 
express the corresponding GDH SRs and their first derivatives. Analogous to the QED case, we 
obtain : 

(i) the GDH SRs: 

oo oo 

(ii) the linearized SRs (valid to leading order in the coupling g): 

oo oo 

fduj (p) 4vra 1 f du (n) 



oo 

4-7ra 1 [ doj 



M 2 " P 
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(iii) the consistency conditions (valid to leading order in the coupling g): 



\ 0=i/^A.<;>. (23c) 



The first derivatives of the cross-sections that enter in Eq. I|23|). to leading order in g, arise through 
the interference of Born graphs Fig. QJa) with the graphs in Fig. H,b) and we find: 

Acrg? = Aa^° p) + Aa% +n) = {2*A[4 + (1 - 2a) (2 + s + 2s)] + 2s\(a + 2) 

- Aln^ 1 4 + (2as + l-x)ln^4l, (24a) 
p — A a — A J 

= + = ^{2A(2 + 2,-,-) + ^ln^-ln^}, (24b) 

Aag.A^ + A.Sf"' = *g,{ ta £±* + »(x/» -*,)}, (24c) 

Aa£> S A„<;""> + Aa<;-" = |^j(2 s - a -x)ln^ + 2A(*-25)}. (24d) 

Using the latter two expressions we easily verify the consistency conditions given in Eq. (|23c|) . 
while, employing the linearized SRs, we obtain: 



oc 



Tre 2 J u lp 



g 2 / /i (4-ll M 2 +3^ 4 ) 2 2 , 3 v 1 

(4^)2 j 1 ~/f _ 1^2 arccos--6^ + 2fM (-5 + 3/T) ln/xj , (25a) 



ooop) = ^ /^ A „<») = rV I , _ M2-/. 2 ) _ t 



= ^y^ A <»=(# 2-«p-ccos f -2^^ . (25b) 

We have checked that Eq. (|25|) agrees with the one-loop calculation done by using the standard 
Feynman-parameter technique. It is interesting that, to this order, the pseudoscalar pion-nucleon 
coupling gives exactly the same result. 

On the other hand, this result does not agree with the covariant ChPT calculation of Ref. [6], which 
are based upon the "infrared-regularization" procedure of Becher and Leutwyler. The discrepancy is 
apparently due to the fact that the "infrared-regularized" loop amplitudes do not satisfy the usual 
dispersion relations. Their analytic properties in the energy plane are complicated by an additional 
cut due to explicit dependence on y/s. In other words, they do not obey the analyticity constraint 
which is imposed on the sum rule calculation. 

It is instructive to examine the chiral behavior of the one-loop result for the nucleon magnetic 
moment. Expanding Eq. (|25|) around the chiral limit (m n = 0), which incidentally corresponds here 
with the heavy-baryon expansion, we have 

4 l0 ° P) = ^{l-27r M -2(2 + 51n^)/, 2 + ^ / i 3 + 0( / i 4 )}, (26) 
4 lo ° p) = j£y {-4 + 2^-2(l-21n^ 2 -^ 3 + 0(/)}. (27) 
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The term linear in pion mass (recall that fi = m n /M) is the well-known leading nonanalytic (LNA) 
correction. On the other hand, expanding the same expressions around the large m n limit we find 



4 loop) = jhf^-^^ +0 ^ ( 28 ) 

4 lo ° p) = tJ^I 2(3-41n/,)-i + 0^- 4 ). (29) 

What is intriguing here is that the one-loop correction to the nucleon a.m.m. for heavy quarks behaves 
as l/rriquark (where m quar k ~ m^.), precisely as expected from a constituent quark-model picture. Here 
this is a result of subtle cancellations in Eq. (|25|) taking place for large values of m n . In contrast, 
the infrared regularization procedure [6] gives the result which exhibits pathological behavior with 
increasing pion mass and diverges for = 2M. 

Since the expressions in Eq. Q25JI have the correct large m n behavior they should be better suited 
for the chiral extrapolations of the lattice results than the usual heavy-baryon expansions or the 
"infrared-regularized" relativistic theory. This point is clearly demonstrated by Fig. |5J where we plot 
the m^-dependence of the full [Eq. (|25|)]. heavy-baryon, and infrared-regularization [6] leading order 
result for the magnetic moment of the proton and the neutron, in comparison to recent lattice data [7]. 
In presenting these results we have added a constant shift (counter-term kq) to the magnetic moment, 
i.e., 

Hp = (l + K 0p + 4 loop) )(e/2Af), (30) 
fin = ( K0 n + 4 lo ° p) )(e/2M) (31) 

and fitted it to the known experimental value of the magnetic moment at the physical pion mass, 
Hp ~ 2.793 and [i n ~ —1.913, shown by the open diamonds in the figure. For the value of the ttNN 
coupling constant we have used g 2 /An = 13.5. The m^-dependence away off the physical point is then 
a prediction of the theory. The figure clearly shows that the SR results, shown by the dotted lines, is 
in a better agreement with the behavior obtained in lattice gauge simulations. 

It is therefore tempting to use the SR results for the parametrization of lattice data. For example, 
we consider the following two-parameter form: 

= {tt^ + '^%- (32b) 

where kq p and Ron are fixed to reproduce the experimental magnetic moments at the physical m n . 
The parameter a can be fitted to lattice data. The solid curves in Fig. |2] represent the result of such 
a single parameter fit to the lattice data of Ref. [7] for the proton and neutron respectively, where 
Op = 1.6 /M 2 and a n = 1.05/M 2 , M is the physical nucleon mass. 
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Figure 2: Chiral behavior of proton and neutron magnetic moments (in nucleon magnetons) to one 
loop compared with lattice data. "SR" (dotted lines): our one-loop relativistic result, "IR" (blue 
long-dashed lines): infrared-regularized relativistic result, "HB" (green dashed lines): LNA term in 
the heavy-baryon expansion. Red solid lines: single-parameter fit based on our SR result. Data points 
are results of lattice simulations. The open diamonds represent the experimental values at the physical 
pion mass. 



3 The nucleon polarizability puzzle 

It is well-known that the Heavy-Baryon ChPT (HBChPT) at order p 3 gives a remarkable prediction 
for the electric and magnetic polarizabilities of the nucleon: 

where gA — 1-26, ~ 93 MeV (related to the ttNN coupling constant used in the previous section 
via the Goldberger-Triemann relation: gA/ fir = g/M). Remarkable, because this is a true prediction 
of HBChPT (there are no counter-terms at this order) and because it appears to be in a very good 
agreement with experiment. 

For now, I am concerned only with the sum of these polarizabilities. On the experimental side, a 
recent determination from the Baldin's sum gives [8]: 

(a E + Pm) p = (13.69 ± 0.14) x 1(T 4 fm 3 , 

(cue + 0M)n = (14.40 ± 0.66) x 10~ 4 fm 3 , (34) 

for proton and neutron, respectively. 

It is also well-known that the A(1232)-resonance excitation gives a large effect to the magnetic 
polarizability. To quantify this effect we use the following Lagrangian for the 7./VA coupling [9, 10]: 

= 2M( m + M^ T * ( i9Mp " U ~ 9E ^ F ") 9 ^ + HX " (35) 
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where \I/ is the nucleoli field, ip^ is the isospin-3/2 spin-3/2 vector-spinor field of the A-isobar, T3 is the 
isospin iVA transition matrix. The coupling constants can be deduced from the empirical knowledge 
of the 7./V — ► A transition strength. Based on the Particle Data Group values for Ml and E2 we 
estimate gu — 3 and gE — — 1- Computing the sum of the s- and ti-channel A contributions, Fig- El 





Figure 3: The A-excitation graphs. 



to the polarizabilities one finds (see [10] for more details): 

(A) 2aff| 







(m + M A y 



0.1 x 10~ 4 fm 3 



/3 (A) 



2ag 2 



M 



1 



7.3 x 10" 4 fm d 



(36a) 
(36b) 



(M + M A ) 2 M a - M 

So, while the HBChPT without Delta's is in a good agreement with experiment, HBChPT with 
Delta's is not at all. Puzzling... One would expect that including the A would improve the situation, 
extend the limit of applicability of our theory higher in energy, into the resonance region. 

There are suggestions voiced now and then that the effect of the A is canceled by the cr-meson 
exchange, or correlated two-pion exchange. In EFT language this corresponds to canceling the A 
contribution by an effect which formally is of higher order in power counting than the A contribution. 
This kind of scenario has recently been explored in Ref. [11] where counter-terms of 0(p 4 ) were 
"promoted" to lower order in order to cancel ffffl . 

Here I would like to argue that possibly there is a more natural explanation within the relativistic 
chiral EFT. To find the leading order prediction of chiral loops relativistically we have computed the 
unpolarized total cross-sections, corresponding to the Born graphs of single-pion photoproduction [12]: 





ttC 
Mlo 3 




2ttC 
Mlo 3 




= 0, 




2vrC 

Mlo 3 



l [u 2 - fi 2 a s] In + 2A [io 2 (a - 2)1 + s^ 2 } X 

+ 2A (aco 2 + s/i 2 



-pSfl In ■ 



uj 2 In 



(3- A 
a + A 



fi 2 (sf3- n z M 2 ) In 



2^2x , /3 + A a + A 



a — A (3 — A a — A 

Substituting these expressions into the Baldin SR, Eq. (jSJ), we obtain: 

406 - 737/x 2 + 304/ - 36/u 6 



(37a) 
(37b) 



isxe + p M i RLo) 



eV 



[3(1 - 4/i 2 + 2/i 4 ) + ±/i 2 ] ln/i + 



+ 



16tt 3 M 3 1 

44 - 788/i 2 + 1500// 4 - 899/i 6 + 215/ - 18^ 10 



6(4 - fi 2 ) 2 



(a E + Pm) 



(RLO) 



e 2 g 2 
16vr 3 Af 3 



In fi + 



3/i(4 - / u 2 ) 5 / 2 
1 



(38a) 



(4 - /i 2 ) 2 



2(2-3, 2 )(ll-5, 2 )-3^ [J^l + a 

3 At ^/4^ V 



(38b) 
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Note that the same result is obtained in the conventional one-loop calculation [13]. The semi- 
relativistic (or, in this case also, chiral) expansion goes as follows: 

n ^(R^n^ e 2 o 2 H ( 48(4 + 31nu) 1521 „ \ . , 



„ , (m n\ eV 11 / 4(1 + 12 In it) H7 ? 



(4vr) 2 M 3 48/x \ llvr 
eV 11 / 4(1 + 121 
(4vr) 2 M 3 48^ \ + liT 

or, numerically (using g 2 /4ir = 13.8, M = 0.9383 GeV, fi = 0.148), 

(qe + Pm) ( p RLO) = 14.5 -5.2 -5.5 + ... = 5.3 (39c) 
{a E + j3 M ) {RLO) = 14.5 -5.5- 0.4 + ... = 8.7 (39d) 

in the usual units. (The total values are consistent with L'vov's numerical calculation [14], if we use 
<7 2 /4-7r = 14.2. Prom that calculation it is clear that a lot of the reduction in the value of the sum 
affects the magnetic polarizability.) 

Therefore, as one can see, the fully relativistic leading order result is substantially different from the 
non-relativistic (heavy-baryon) limit. The relativistic corrections which are suppressed by m^/M ~ 
1/7, and hence are supposed to be small, appear with large coefficients and actually are not that 
small. The good news is that this apparently allows us to accommodate the large effect of the A 
isobar. In fact, the A contribution now improves the agreement with experiment. Adding the RLO 
and A numbers we have: 

(a E + Pm) ( p RLO+A) = 5.3 + 7.2 = 12.5 x 10" 4 fm 3 (40) 
(a E + Pm) {RLO+A) = 8.7 + 7.2 = 15.9 x 10~ 4 fm 3 . (41) 

At this order there is also an effect of ttA chiral loops, but those are not yet computed relativistically. 



4 Conclusion 

The chiral EFT of QCD provides a description of the low-energy hadronic reactions and that allows one 
to extract hadron properties from experiments. On the other hand, it predicts the chiral behavior of 
these properties and that allows one to make a link to the lattice QCD calculations. These are the two 
fronts which at present make the chiral EFT indispensable in relating QCD to low-energy observables. 
The purpose of this talk is to demonstrate, on simple examples of nucleon magnetic moment and 
polarizabilities, that manifestly relativistic calculations do a better job than the "heavy-baryon" ones 
on both fronts. 

The calculations done in this work were based on the real-Compton-scattering sum rules, such as 
GDH and Baldin sum rules. However, the results are not different from what one would obtained in 
the usual loop calculations, provided no manipulations (e.g., infrared regularization) which change the 
analytic structure are made. As is shown in the works of Gegelia et al. [15], there is no problem with 
power-counting in this, straightforward, formulation of covariant ChPT, if the renormalization scale 
is set in a suitable way. 



10 



Acknowledgements 

I would like to extend my gratitude to the organizers for the invitation, to all the younger crowd for 
the great time, and of course to the Sicilian mob for sparing our lives despite some tensions over our 
dining preferences. 

This work is supported in part by DOE grant no. DE-FG02-04ER41302 and contract DE-AC05- 
84ER-40150 under which the Southeastern Universities Research Association (SURA) operates the 
Thomas Jefferson National Accelerator Facility. 



References 

D. Drechsel, B. Pasquini and M. Vanderhaeghen, Phys. Rept. 378 (2003) 99. 



D. Drechsel and L. Tiator, arXiv:nucl-th/0406059 ; see also contributions of R. Beck and P. Grab- 
mayr in these proceedings. 

G. Altarelli, N. Cabibbo, and L. Maiani, Phys. Lett. B 40 (1972) 415. 
D.A. Dicus and R. Vega, Phys. Lett. B 501 (2001) 44. 

V. Pascalutsa, B. R. Holstein and M. Vanderhaeghen, Phys. Lett. B 600 (2004) 239. 
B. Kubis and U. G. Meifiner, Nucl. Phys. A 679 (2001) 698. 

J. M. Zanotti, S. Boinepalli, D. B. Leinweber, A. G. Williams and J. B. Zhang, Nucl. Phys. Proc. 
Suppl. 128 (2004) 233 |arXiv:hep-lat/0401029| . 

D. Babusci, G. Giordano and G. Matone, Phys. Rev. C 57, 291 (1998). 

V. Pascalutsa, Phys. Rev. D 58 (1998) 096002; V. Pascalutsa and R. Timmermans, Phys. Rev. C 
60 (1999) 042201. 

V. Pascalutsa and D. R. Phillips, Phys. Rev. C 67 (2003) 055202; Phys. Rev. C 68 (2003) 055205; 
|arXiv:nucl-th/0308065 



R. P. Hildebrandt, H. W. Griesshammer, T. R. Hemmert and B. Pasquini, Eur. Phys. J. A 20 
(2004) 293; see also contribution of H. W. Griesshammer in these proceedings. 

B. R. Holstein, V. Pascalutsa and M. Vanderhaeghen, in preparation. 

V. Bernard, N. Kaiser and U.-G. Meifiner, Nucl. Phys. B 373 (1992) 346; A. Metz and D. Drechsel, 

Z. Phys. A 356 (1996) 351. 

A. I. L'vov, Phys. Lett. B 304 (1993) 29. 

J. Gegelia, G. Japaridze and X. Q. Wang, J. Phys. G 29 (2003) 2303 |arXiv:hep-ph/9910260| ; 
T. Fuchs, J. Gegelia, G. Japaridze and S. Scherer, Phys. Rev. D 68 (2003) 056005. 



11 



